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' Abstract We present spectroscopic observations of galaxies in 4 clusters at z = 0.7 - 0.8 and in one cluster at z ~ 0.5 obtained 

with the FORS2 spectrograph on the VLT as part of the ESO Distant Cluster Survey (EDisCS), a photometric and spectroscopic 
survey of 20 intermediate to high redshift clusters. We describe our target selection, mask design, observation and data reduction 
' procedures, using these first 5 clusters to demonstrate how our strategies maximise the number of cluster members for which 

I we obtain spectroscopy. We present catalogues containing positions, I-band magnitudes and spectroscopic redshifts for galaxies 

in the fields of our 5 clusters. These contain 236 cluster members, with the number of members per cluster ranging from 30 
to 67. Our spectroscopic success rate, i.e. the fraction of spectroscopic targets which are cluster members, averages 50% and 
ranges from 30% to 75%. We use a robust bi weight estimator to measure cluster velocity dispersions from our spectroscopic 
redshift samples. We also make a first assessment of substructure within our clusters. The velocity dispersions range from 
400 to 1100 kms '. Some of the redshift distributions are significantly non-Gaussian and we find evidence for significant 
substructure in two clusters, one at z ~ 0.79 and the other at z ~ 0.54. Both have velocity dispersions exceeding 1000 kms"' 
but are clearly not fully virialised; their velocity dispersions may thus be a poor indicator of their masses. The properties of 
these first 5 EDisCS clusters span a wide range in redshift, velocity dispersion, richness and substructure, but are representative 
of the sample as a whole. Spectroscopy for the full dataset will allow a comprehensive study of galaxy evolution as a function 
of cluster environment and redshift. 
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* Based on observations obtained at the ESO Very Large Telescope 
(VLT) as part of the Large Programme I66.A-0I62 (the ESO Distant 
Cluster Survey). 

** Table |4| is only available in electronic form at the CDS 
via anoymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
^http://cdsweb.u-strasbg.fr/ Abstract.html 



1. Introduction 

Galaxy clusters are convenient and popular sites for study- 
ing environmental influences on galaxy evolution and the 
systematic analysis of large cluster samples is an important 
step towards a quantitative understanding of the processes 
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involved. Large area cluster searches began with the pio- 
neering pho tographic s u rveys of lAbelll lll958h . I^vickv et al] 
ill 9681) and LShectmanI lll985ll. Mor e recently, surve ys for 
high-redshif t clusters (e.g^Gun n et al.>il98 6: Gioia et al. 19901 
Couch etalJlT9 91: Postman et al .ll996H Gladders & Yee 2000; 
Gonzalez et al.l2001l) have begun to provide large enough sam- 
ples to allow an exploration of the evolution of clusters and the 
galaxies within them out to significant lookback times. 

Spectroscopic surveys of the full galaxy population of dis- 
tant galaxy clusters are few in number. At intermediate red- 
shift they include the studies of Couch & Sharpies ( 1987), and 
of the CNOC and MORPHs collaborations (Yeeetal. 1996; 
iBaloghet all 19971; .Dressier et al.,1999.; Poggianti et aL 1 999 ), 
as well as the studies of individual clusters by van Do kkum, 
Kelson, Tran and collaborators (e.g. 'van Dokkum et al Booot 
iTran et allEool) . At high er redshift there is the survey by 
Postrn an and collaborators dPostman et alJl998llPostman et alJ 
1200 ih and the study of the single very rich cluster MS 1054- 
03 bv lvan Dokkum et"al..i.l999.. Such surveys have contributed 
substantially to our understanding of the star formation histo- 
ries of cluster galaxies, revealing strong evolution in their av- 
erage star formation activity since z ~ 0.5, thus over the last 5 
Gyr At redshifts beyond 0.5, too few clusters have been stud- 
ied so far to draw strong conclusions about galaxy evolution 
and the dependence of galaxy properties on environment (e.g. 
van Dokkum et al. 1999, 2000; Postman et al. 1998, 2001). 

Measurements of cluster velocity dispersions are ofte n used 
to estimate cluster masses (e.g. iFisher et al. . .,1998.: .Tran et alJ 
Il999t iBorgani et al.lll999t iLubin etaLll2002l). Early mea sure- 
ments of velocity dispersions (e.g. Danese et al. 1980) usu- 
ally assumed a Gaussian velocity distribution. As datasets 
have increased in quality and size, significant deviations from 
Gaussia n behaviou r have often been found (e.g. Zabludof f et all 
Il993.; ■Fisher et alJil998). These are likely associated with the 
substantial structural irregularities (usually termed substruc- 
ture) which are observ ed in many clusters (Geller & Beers 
1982: IPressler & Shectman 1988,). Any statistic used to es- 
timate cluster velocity dispersion should therefore be robust 
against outliers and against vari ations in t he shape of the un- 
derlying velocity distribution (Beers et al.. 1990) . Cluster ve- 
locity dispersion measurements have been p ubhshed for only 
a small n umber of high redshift clusters (e.g. Tran et al. 1999; 
iLubin et aL 2002) and analyses of cluster substructure at such 
high redshift are scarce. In this paper we present velocity dis- 
persions and substructure analyses for 5 of the high redshift 
clusters in the ESO Distant Cluster Survey (EDisCS). 

The ESO Distant Cluster Survey is a photometric and 
spectroscopic survey of 20 clusters selected from th e Las 
Campanas Distant Cluster Survey (LCDCS) ( Gonzalez et alJ 
I2OOII) . Clusters in the LCDCS were selected to be amongst the 
most significant surface brightness peaks in a smoothed drift 
scan map of 1 10 square degrees of the southern sky. One half 
of our sample was chosen to have high estimated redshift in 
the LCDCS {Zesi ~ 0.8) while the other half was chosen to 
have intermediate estimated redshift (zesr ~ 0.5). The appar- 
ent reality of these candidates was established before inclusion 
in the final EDisCS sample by moderate depth VLT imaging 
in two passbands. This allowed verification of the expected 



overabundance, colour and magnitude of red sequence galax- 
ies (lOonzalez et al. 2002). 

Our EDisCS observing programme is now complete. It in- 
volves deep optical photometry using FORS2 at the VLT (14 
nights), near-IR photometry with the SOFI instrument at the 
NTT (20 nights), and multi-slit spectroscopy using the FORS2 
spectrograph at the VLT (22 nights). Our optical photometry 
comprises V, R and I imaging for the 10 high-redshift clusters 
and B, V and I imaging for the 10 intermediate-redshift clus- 
ters (White et al. 2004, in preparation). The field of view of the 
F0RS2 instrument is 6.8' x 6.8'. After dithering, the field of 
view with the maximum depth of exposure in our imaging data 
was approximately 6.5' x 6.5'. Near-IR data in J and K will 
be presented in Aragon-Salamancaet al. 2004 (in preparation). 
The field of view of our SOFI images, taking account of the 
dithering and overlapping of exposures, was 6.0' x 4.2' for our 
clusters at z~0.5 and 5.4' x 4.2' for our clusters at z~0.8: see 
Aragon-Salamanca et al. 2004, in preparation, for details. This 
ground-based dataset has been complemented in February-June 
2003 by 80 orbits of HST ACS imaging for our 10 high-redshift 
clusters (Desai et al. 2004, in preparation). 

The EDisCS project aims to investigate galaxy morphol- 
ogy, structure, masses, stellar populations and star formation 
histories as a function of both redshift and cluster environment, 
and to study cluster mass and structure. Our spectroscopy al- 
lows cluster membership to be defined using precise redshifts 
and it constrains the evolutionary histories of galaxies in dense 
environments by providing estimates of (or at least constraints 
on) star formation rates and histories, chemical abundances, in- 
ternal kinematics, and stellar and dynamical masses. We are 
studying both early- and late-type galaxy evolution by mea- 
surement of the Fundamental Plane and Tully-Fisher scaling re- 
lations as a function of redshift and environment. The ages and 
metallicities of cluster absorption-line galaxies are constrained 
using absorption line-strength index measurements and evolu- 
tionary population synthesis models. The star formation histo- 
ries of cluster galaxies are being studied using absorption line 
and emission line information, as well as broad-band spectral 
energy distributions and HST/ACS morphologies. 

In this paper we present the procedures we have adopted 
to carry out the spectroscopic part of the EDisCS programme, 
and we give first results for a set of five clusters. The main fo- 
cus of this paper is a presentation of the spectroscopic data for 
these clusters, one of which is at z ~ 0.5 (CI 1232.5-1250) 
while the other four are at z = 0.7 - 0.8 (CI 1040.7-1155, 
CI 1054.4-1146, CI 1054.7-1245 and CI 1216.8-1201). We 
will abbreviate these names as CI 1232, CI 1040, CI 1054-11, 
CI 1054-12 and CI 1216, respectively, in the following. In 
Section 12 we describe our target selection strategy, our mask 
design and our VLT FORS2 observations. In Section |3] we 
describe the reduction of these VLT FORS2 observations. In 
Section |4] we present galaxy redshift measurements and spec- 
troscopic catalogues, and we assess our observational com- 
pleteness and the success of our observing strategy. In Section 
|5]we present velocity dispersion estimates for our clusters. In 
Section |6l we study cluster substructure using maps of galaxy 
positions and also the Dressler-Shectman (1988) method. A 
summary is provided in Section^ 
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Table 1. The EDisCS spectroscopic observing runs 



Run 


Dates 


Number of nights 




(noon-to-noon) 


Total 


Usable 


Run 1 


06 Feb 2002 - 09 Feb 2002 


3 


1.6 


Run 2 


11 Mar 2002- 19 Mar 2002 


8 


6.9 


Run 3 


26 Mar 2003 - 06 Apr 2003 


11 


11.0 


Run 4 


23 Mar 2004 - 26 Mar 2004 


3 


3.0 



2. Observations 

Spectroscopic observations were completed during the four ob- 
serving runs listed in Table [2 During our first run we took rel- 
atively short exposures to establish the reality of our clusters, 
to measure their redshifts and to tune our photometric redshifts 
to allow stricter selection criteria for subsequent runs. During 
our second run deep spectroscopic exposures were completed 
for the 5 EDisCS clusters discussed in this paper, 4 high red- 
shift clusters and 1 intermediate redshift cluster. Short expo- 
sures were acquired in addition in this second run for EDisCS 
clusters for which data were not obtained in the first run due to 
bad weather. During the third run, one or two additional masks 
were observed for the 5 clusters presented in this paper, and 
long-exposure data were obtained for most of the other EDisCS 
clusters. Finally, in the fourth run we completed our spectro- 
scopic programme for the full sample, obtaining data for clus- 
ters not discussed in this paper. 

In section IZTI we describe the target selection strategy. In 
section IT^ we describe the VLT FORS2 MXU mask prepara- 
tion. In section l23l we describe the observations themselves. 

2.1. Target Selection Strategy 

When selecting the spectroscopic targets, the main aim was 
to produce a well-defined sample of cluster galaxies that was 
as statistically complete as possible. Given the available imag- 
ing data, selection in the I-band was the obvious choice: such 
red-band selection reduces biases towards star-forming or un- 
extincted galaxies. Our parent sample consisted of all objects 
with 18.6 < / < 22 for the z ~ 0.5 cluster candidates and 
19.5 < / < 23 for the z ~ 0.8 ones. Since we are com- 
pleting spectroscopy using 1" slits, for our target selection 
we used 1" radius aperture magnitudes. The bright magni- 
tude cut was imposed to reject galaxies that are too bright 
to be cluster members at the expected redshifts of our clus- 
ters. These bright magnitude limits are ~ 1 mag brighter than 
brightest cluster galaxies at z ~ 0.5 and z ~ 0.8 respectively 
jAragon-Salamanca et al.lfl998h . and are thus very conserva- 
tive. 

Such straightforward I-band selection produces a sample 
where the selection criteria are easy to model and to reproduce, 
but where the contamination by non-cluster galaxies and by 
stars is quite high. For example, for the extremely rich z = 0.83 
cluster MS 1054-03, Keck spectroscopy for an I-band selected 
sample found that only ^ 40% of the gala xies with measured 
redsh ifts were actually cluster members ( Ivan Dokkum et all 
I2OOO ). Our optically-selected clusters have a range of rich- 
nesses, so we would expect our success rate for a purely I- 



band selected sample to be substantially worse than this in 
most cases. To avoid wasting telescope time, we decided to use 
the available imaging information to impose additional selec- 
tion criteria. The aim is to reject a significant fraction of non- 
members while keeping the sample of cluster galaxies as close 
as possible to an I-band selected one. Since completeness of 
our cluster samples is our main concern, our rejection criteria 
were chosen to be conservative. 

Our target rejection criteria were based primarily on pho- 
tometric redshifts obtained with a modified version of the code 
hyperz' (Bolzonella et al. 2000). Photometric redshifts have 
proven to be quite reliable when a reasonably-complete set of 
photometry, including near-infrared bands, is available. When 
tested against spectroscopic reds hifts, the typical errors in the 
photometric redshifts are ^ 0.1 ( Bolzonella et al . 2000; Pello 
et al. 2004, in preparation). The hyperz code is based on a 
standard SED fitting minimization procedure, broad-band fil- 
ters allowing the detection of strong spectral features, such as 
the 4000 A break, the Lyman break or strong emission lines. 
The template library used includes 6 evolutionary synthetic 
SEDs computed using the Bruzual & Chariot code (Bruzual 
& Chariot 1993), matching the observed colours of local field 
galaxies at z~ 0, from E to Im type, in addition to a set of em- 
pirical SEDs compiled by Coleman, Wu and Weedman (1980), 
and 2 starburst galaxies from the Kinney et al. (1996) template 
library (SBl and SB2). The internal reddening law is taken 
from Calzetti (2000), and considered as a free parameter with 
Ay ranging between and 1.5 magnitudes (E(B-V) between 
and ~ 0.37 mags). For each cluster field a correction for galac- 
tic extinction was included using the value of E(B-V) derived 
from Schlegel et al. (1998) for the cluster centre. Galactic E(B- 
V) corrections in the EDisCS fields typically range between 
0.03 and 0.08 magnitudes. The reader is referred to the refer- 
ence paper (Pello et al. 2004, in preparation) for a more detailed 
description. 

Photometric redshift accuracy depends on the wavelength 
domain covered by the photometric SEDs, i.e. the filter set. 
The set of filters used by EDisCS was designed to cover the 
relevant wavelength domain in the restframe of mid z~0.5 and 
high z~0.8 clusters, in particular to bracket the 4000 A break 
within the relevant redshift intervals. As explained in White et 
al. (2004), deep optical photometry was obtained using FORS2 
at the VLT, in BVI and VRI bands respectively for our mid-z 
and the high-z clusters. In addition, deep near-IR images were 
obtained with SOFI at NTT, in K' and JK' bands respectively 
for our mid-z and the high-z clusters (see Aragon-Salamanca et 
al. 2004, in preparation for details). These optical and near-IR 
band data were used to derive our photometric redshifts. 

After imposing the magnitude cuts described above, targets 
were selected if their photometric redshift obeyed Zcius - 0.2 < 
Zphot ^ Zcius + 0.2, where Zcius was the estimated cluster red- 
shift. The cluster redshift was estimated, initially, from the 
peak in the photometric redshift histograms (see Pello et al. 
2004, in preparation). After our first short-exposure masks 
were observed and analysed, cluster redshifts based on our 
spectroscopy were used instead. 

' http://webast.obs-mip.fr/hyperz/ 



4 



C. Halliday et al.: Spectroscopy of clusters in the ESO Distant Cluster Survey (EDisCS) 



As an additional measure to ensure that no likely cluster 
members were dropped from the sample, we took into account 
that the photometric redshift code provides not only a "best 
guess" photometric redshift, but also a redshift probability dis- 
tribution for each object. Objects with Zphot outside the Zcius±0.2 
interval were selected for spectroscopy if the ;i'^-based proba- 
bility of the best fit provided by the photometric redshift code 
at the cluster redshift was larger than 50%. Because the red- 
shift probability distributions were reasonably peaked for most 
objects, this only increased the number of selected targets by 2- 
3%. For this reason, the exact value of the probability threshold 
for selection had very little eff'ect on our target list (see Pello et 
al. 2004, in preparation for details). 

Finally, the photometric redshift code also provides some 
information to help with star-galaxy separation. For each ob- 
ject, the code computes two independent values, A^g and Nt, 
based on how well the photometry corresponds to a galactic or 
a stellar SED (Pello et al. 2004). A^g = means "the object 
is almost certainly not a galaxy"; A^g - 1 "the object could 
be a galaxy"; and A^g = 2 "the object is almost certainly a 
galaxy". Similarly, A^, can have values 0, 1 and 2 with similar 
meanings but for stars. Objects with A^g = and > I were 
excluded from the spectroscopic sample as very probable stars. 
As a safety precaution, we have checked that the vast major- 
ity (at least ^ 95%) of the objects rejected as probable stars 
based on their photometry have stellar PSFs, while virtually 
all the rest are blended objects. Moreover, each spectroscopic 
mask had slits placed on several stars (typically 3-5) for ac- 
quisition purposes. These stars were selected from the objects 
photometrically classified as very probable stars. There were 
52 such objects in the masks of the five clusters discussed in 
this paper, and the spectroscopy confirmed all of them as stars. 

Indeed, the photometric criteria used to remove stars from 
the target list is more conservative than using the PSF alone: a 
significant number of objects with stellar PSFs remained in the 
target list in the 2002 observing runs, and they were confirmed 
as stars by the spectroscopy. To further reduce our stellar con- 
tamination, in the 2003 and 2004 observing runs we decided to 
reject in addition objects which have clear stellar PSFs in the 
I-band images. This decision was based on the excellent image 
quality of our I-band data (median seeing 0.6", cf. White et al. 
2004), and on the experience gained in the 2002 runs. For each 
cluster, a set of 20-30 bona-fide stars was selected, their ellip- 
ticities and FWHMs measured, and the mean ((FWHM)) and 
standard deviation of the FWHMs ((Tpwhm) determined. Then, 
objects were removed from the target list of that cluster if their 
FWHM was less than (FWHM)-i-2o-pwHM and their ellipticities 
were less than 0. 1 . Typically, (FWHM) ^ 0.6" and ctfwhm - 
0.05". This procedure removed only ~ 5% of the candidate 
targets. Note that the results presented in this paper are based 
mostly on spectroscopy obtained in 2002, and thus this addi- 
tional selection criterion was not applied there. Had we applied 
it, out of the 231 spectroscopically-confirmed cluster mem- 
bers, not a single one would have been rejected on the basis 
of its PSF. Indeed, not a single object with spectroscopically- 
confirmed redshift larger than ^ 0.0 would have been rejected 
on the basis of its PSF. 



In addition to these photometric criteria, the mask design 
imposes geometric constraints in order to avoid overlapping 
spectra and to ensure the necessary wavelength coverage (see 
next section). In the mask design process, it was sometimes 
impossible to find targets that fulfilled all the selection crite- 
ria at a given mask location. In such cases, in order to fill the 
masks, slits were placed on an ad hoc basis on objects that had 
not been selected as cluster galaxy candidates, and therefore 
were expected to be field galaxies based on their photometric 
redshifts. A test of how well our selection worked is provided 
by the number of such photometrically-selected field galaxies 
that turned out to be cluster members. For the five clusters dis- 
cussed in this paper, redshifts were obtained for 65 such galax- 
ies, and none of them turned out to be a cluster member Thus 
we are reasonably confident that our selection criteria have ex- 
cluded, at most, 1-2% of real cluster galaxies above our I-band 
limit. Once the spectroscopy for our complete EDisCS data set 
is analysed, we will have tighter constraints on this number. 

A measure of the effectiveness of our selection procedure is 
that our success rate, measured as the fraction of spectroscopic 
targets that turned out to be cluster members, is, on average, 
- 50% for the five clusters described in this paper. If we had 
used a pure I-band selection, the success rate would have been 
under 30%, and our use of 8-m telescope time substantially 
less efficient. In summary, we are quite confident that our target 
selection process yielded a cluster galaxy sample that is very 
close to a purely I-band selected one, but with almost twice as 
many spectroscopically-confirmed cluster members as would 
have been the case if we had not made use of photo-z rejection 
of non-members. 

2.2. Automatic Creation of the Masks 

We developed a program to automatically place MXU mask 
slits on targeted galaxies (iPomer.2004.) . 

Mask design was performed using our EDisCS I-band im- 
ages, since these best correspond to the wavelength domain 
of the grism 600RI+19 which we had chosen for our spec- 
troscopy. We fixed the right ascension axis to be parallel to the 
X axis of the detectors. The centre of the mask was chosen in 
order to maximize the number of galaxies to be observed. For 
the 2002 observations we matched the spectroscopic field cen- 
tre with the position of the brightest cluster galaxy (BCG) in X 
and the center of the photometric image in Y. For 2003, we had 
to take into account the change of F0RS2 CCDs and in partic- 
ular a 5"5 gap between the two new detectors. We allowed a 10 
to 20" shift along the Y axis. 

Reference stars were selected to aid the calculation of shifts 
in RA, Dec and rotation required to precisely align each MXU 
mask on the sky. Seven reference stars were chosen for each 
mask to have 17.5 mag < I < 19.9 mag and to be uniformly 
distributed across the field of view. Short (5') slits were placed 
on 2-3 stars to help optimise the final mask alignment and to 
allow seeing to be measured during observations. 

Our science aims required the observation of a minimum 
spectral region including the [OII]/l3727 and H(J/14101 spec- 
tral fines for the cluster galaxy candidates. This restframe do- 
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main was shifted for each cluster, according to its photometric 
redshift for the preparation of our first February 2002 spectro- 
scopic run, or using the spectroscopic cluster redshift for the 
preparation of our long-exposure masks. 

The slit width was 1". We selected slit lengths of 1" for 
cluster galaxy candidates, 6" for field galaxy candidates, and 
5" for star reference slits. For the cluster BCG galaxy a 10" 
slit was not sufficiently long and we chose a slit length equal 
to lOxFWHM. We fixed a minimum separation between slits 
of 1" to avoid any superposition of spectra. The galaxies were 
centered in the slits when possible, otherwise at least 3" away 
from the slit edges to ensure good sky substraction. 

In 2002, slits on cluster galaxy candidates were aligned 
with the galaxy major axis when hyperz indicated a spiral 
type spectral energy distribution. The position angle provided 
by SExtractor was used as long as it did not exceed 45 degrees. 
A vertical slit was selected for all remaining galaxies. After 
reduction of the 2002 runs, we identified a number of emis- 
sion line galaxies which had been classified as ellipticals by 
hyperz. In subsequent runs slits were aligned with the major 
axis for all galaxies provided that the slit inclination angle did 
not exceed 45 degrees. 

Galaxies were selected by order of priority: 1) BCG or the 
closest possible galaxy; 2) confirmed cluster galaxy, or target 
for which the redshift determination was still ambiguous in 
the short exposure spectra; 3) cluster galaxy candidates, based 
on the criteria presented in Section ITTl 4) star slits; 5) field 
galaxy candidates. For each of the categories 1 to 3, the spec- 
troscopic target selection was completed so that the brightest 
cluster member candidates in the considered region of the mask 
were selected first. 

Once all possible primary targets had been exhausted, we 
filled masks with 6" slits placed preferentially on cluster galaxy 
candidates, with no constraint on the wavelength range. Finally 
a visual inspection of aU masks was performed. 

2.3. Observations 

Spectroscopic data were acquired during four observing runs 
using the MXU multi-object mask facility of the FORS2 spec- 
trograph mounted on the VLT Yepun UT4 telescope, ESO 
Paranal, as described in Tables 1 and 3. Details of our obser- 
vational setup are provided in Table |2 

The observing log is given in Table |3] For each of the 20 
EDisCS clusters the total exposure time for each mask ob- 
served in the four observing runs is listed. The individual expo- 
sures were typically 1800 s. The masks are labelled in a simple 
way: masks used in run 1 are called mask 1, masks used in 
run 2 are called mask 2, 3 and 4, masks used in run 3 are called 
mask 5, 6, 7 and 8, and masks used in run 4 are called mask 9, 
10 and 11. Since not all clusters were observed in run 1, not all 
clusters have a 'mask 1', et cetera. 

During run 1 relatively short exposures (typically 1 hour) 
were completed for the 10 clusters of our high-z sample and 
for one cluster in our intermediate sample. Time lost to bad 
weather meant that short exposures could not be obtained for 
the remaining mid-z clusters in this run. For the data acquired, 



Table 2. The MXU multi-slit instrumental set-up 



Run 


run 1+2 


run 3 + 4 


Telescope 


VLT-UT4 


Instrument 




FORS2 


Mode 


MXU (i.e 


mask multi-slit) 


Grism 


GRIS.600RI+19 


Filter 




GG435 


Spectral resolution 


6.e 


lAFWHM 


Detector(s) 


Site 


2xMIT 


Detector size 


2048 X 2048 


2 X (2048 X 1024)" 


Pixel Size 


24/jm X 24yum 


30/jm X SO/jm" 


Readout Noise 


541 e- 


3.15 e- 


Conversion factor 


1.91 e-/ADU 


0.70 e-/ADU 


MXU slit width 


I'.'O 


I'.'O 


Spatial scale 


a'20/pixel 


a'25/pixel'' 


Dispersion 


1.32 A/pixel 


1.66A/pixeP 


Wavelength range'' 


5300-8000 A 


5 120-8450 A 



" After binning the CCDs by 2 x 2, which is the default. 
For a slit at the centre of the field of view. 



galaxy redshifts were measured when possible. Masks for our 
run 2 were designed to target established cluster members and 
additional candidate cluster members for which short expo- 
sures were not obtained or a redshift could not be measured 
due to insufficient signal-to-noise. 

In run 2 deep exposures were obtained for the 5 clusters that 
are the main focus of this paper: CI 1232 from the mid-z sample 
and CI 1040, CI 1054-11, CI 1054-12 and CI 1216 from the 
high-z sample. The typical exposure time per mask was 2 hours 
at mid-z and 4 hours at high-z. In addition short exposures were 
obtained for a further 8 clusters from the mid-z sample. Seeing 
was in general subarcsecond and typically varied between 0'.'7 
and 1'.'7. 

In run 3 long exposures were obtained for an additional 1- 
2 masks for the 5 clusters that aheady had long masks from 
run 2. This completed the spectroscopic data for these clusters. 
During run 3 long exposures were obtained for 8 additional 
clusters, with 2^ masks per cluster totaling 27 masks. In run 4 
the planned long exposures on the high-z sample were com- 
pleted, and semi-long exposures (I hour per mask) were car- 
ried out for clusters from the mid-z sample. These data will be 
presented in detail in forthcoming papers. 

The calibration frames obtained in the 4 runs were as fol- 
lows. During the bright part of twilight sky flat-field exposures 
were acquired. During the darker part of twilight calibration 
stars were observed. Hot stars (spectral type late O to early 
B) for the telluric absorption correction were observed with a 
longslit of width I'.'O. Spectrophotometric standard stars for the 
flux calibration were observed using the FORS2 MOS mode. 
For our purpose this provided a short 'longslit' of width 5'.'0 
that could be placed anywhere in the field of view. We com- 
pleted exposures at the extreme left, the centre and the extreme 
right of the field of view, bracketing the range of slit positions 
and hence wavelength ranges of our MXU galaxy slits. 

Bias, screen flat-field and He+Ar+Ne+HgCd comparison 
arc-lamp exposures were taken by an automatic batch proce- 
dure executed after each night of observation. 5 biases, 6 screen 
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Table 3. Observing log for the four EDisCS spectroscopic observing runs 



Cluster 


Run 1 

ml 




_Run 2 






Run 3 






_Run4 




Repeats 


ml 


m3 


m4 


m5 


m6 


m7 


m8 


m9 


mlO 


mil 


Mid-z clusters: 


























CI 1018.8-1211 


1800 




1500 




5400 


5400 


7200 










ml-m5; m3-m6 


CI 1059.2-1253 




1800 






5400 


7200 


7200 










m2-m5 


CI 1119.3-1129 




1800 














3600 


3600 


3600 


m2-m9 


CI 1202.7-1224 




1800 














3600 


3600 


3600 


m2~m9 


CI 1232.5-1250 




7200 


6360 


7200 


7200 
















CI 1238.5-1144 




1800 














1200 






m2-m9 


CI 1301.7-1139 




1800 














3600 


3600 


3600 


m2-m9 


CI 1353.0-1137 




2400 














3600 


3600 


3600 


m2-m9 


CI 1411.1-1148 




1500 






10800 


7200 


9000 










m2~m5 


CI 1420.3-1236 




1800 














3600 


3600 


3600 


m2-m9 


High-z clusters: 


























CI 1037.9-1243 


3600 








7200 


10800 


10800 


10800 








ml-m5 


CI 1040.7-1155* 


3600 


10800 


14400 


14400 


14400 


14400 












ml-m2 


CI 1054.4-1146* 


3600 


10800 


14400 


14400 


14400 














ml-m2 


CI 1054.7-1245* 


3600 


10800 


14400 


14400 


14400 














ml-m2 


CI 1103.7-1245 


3600 








10800 


10800 






14400 


14400 




ml-m5 


CI 1122.9-1136 


3600 
























CI 1138.2-1133 


3600 








3600 


7200 


7200 


7200 








ml-m5 


CI 1216.8-1201* 


3600 


10800 


14400 


14400 


14400 














ml-m2 


CI 1227.9-1138 


3600 








10800 


14400 


14400 


14400 








ml-m5 


CI 1354.2-1230 


1800 








12600 


14100 


15000 


17400 








ml-m5 



Notes — ml, m2, mil denotes mask 1, mask 2, mask 11. The observing run can be deduced from the mask number. A typical 
observing strategy was to take a short exposure on a cluster in one run, followed in a subsequent run by a longer exposure in which only the 
spectroscopically confirmed cluster members from the short mask were repeated. Such mask pairs are listed in the column 'Repeats'. The five 
clusters marked '*' are the main focus of this paper. These are the clusters for which long exposures were obtained in run 2. 



flat-field exposures and a single arc-lamp exposure were ob- 
tained for each mask and set of stellar observations. 



3. Data reduction 

Data reduction was completed primarily using IRAF^. 
Specialised algorith ms developed for the re duction of similar 
FORS2 MXU data llMilvang-JenseiJl2003h were adopted for 
several of the reduction steps. In the reduction of the data from 
run 3, the images from the two CCDs (cf. TableO were reduced 
separately. We describe the data reduction of spectra acquired 
in runs 1 and 2, and for data from run 3 for the five clusters that 
are the principle focus of this paper 

Bias subtraction was completed by subtracting a combined 
bias image. For the run 3 data an overscan level correction 
was made. Following this cosmic ray hits (hereafter 'cosmics') 
were identified. Different approaches were taken depending 
on the number of consecutive exposures obtained for a given 
mask. For the 'long masks' typically 4 or 8 exposures of 30 min 
were available. In thi s case, cosmics were i dentified using a 
method developed bv lMilvang-Jenseiil ( l2003l) described in de- 
tail in AppendixIXl 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Single exposures were treated for cosmics using a me- 
dian filter of dimensions 5 pixels along the dispersion axis 
and 1 pixel along the spatial axis for data acquired in run 1. 
Double exposures acquired in run 1 were treated by comparing 
each frame and taking the minimum value of each correspond- 
ing pixel value as that of the output frame. Data obtained in 
run 2 with 1-3 succes sive exposures were processed using the 
LAcosmic algorithm ( Ivan Dokkuml200ll) : a very conservative 
choice of rejection level was made to prevent removal of sky- 
lines; some small cosmics remained and were treated interac- 
tively using the IRAF task cosmicrays. 

The spectra in the upper part of the field curve like a 'U' 
while the spectra in the lower part of the field curve like an 
upside-down 'U'. The maximum peak to peak effect was 1.2" 
(i.e. 6px for the 2002 data). This geometrical distortion (S- 
distortion) was mapped using the edges of the individual 'spec- 
tra' in the screen flats. Specifically, the screen flats were con- 
volved with a [-1, 1] kernel after which the absolute value was 
taken. This produced an image in which the edges of the spec- 
tra show up as positive features (the convolution is equivalent 
to moving the image downwards by 1 pixel and subtracting it 
from itself.) These features were then used to map the distor- 
tion. All frames (science, flat fields and arcs) were then cor- 
rected for the distortion. This created spectra parallel to the X 
axis. 

For each mask the combined screen flat exposures were 
used to construct two types of flat fields; a pixel-to-pixel flat 
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field and a flat field in which no normalisation was made in the 
wavelength direction. The pixel-to-pixel flat was constructed 
by normalising each row using a 20-piece cubic spline, i.e. by 
removing the wavelength dependence caused by the combined 
effect of the SED of the screen flat lamp and the transmission 
function of the system (mirrors, grism, CCD). The pixel-to- 
pixel flat field was applied to the corresponding science frame. 
In this way the original counts are preserved, allowing the pho- 
ton noise to be readily computed. The second type of flat field 
was constructed as follows. The combined screen flat image 
for the given mask was cut up and each slit spectrum wave- 
length calibrated. The slits in the masks would be at differ- 
ent jc-positions, and the slit closest to the centre was located. 
Only the overall level in the corresponding screen flat spectrum 
was normalised, leaving the wavelength dependence intact. The 
normalisation constant was chosen to give a level of one at a 
reference wavelength of 6780 A, which is central wavelength 
of the grism according to the manual. This normalisation con- 
stant was used to normalise all other screen flat spectra for the 
mask. These flat field spectra were used to create a second set 
of science frames required for the flux-calibration of science 
spectra (cf. below). 

Since the frames had been geometrically rectified, the in- 
dividual slitlet spectra within each CCD frame were identified 
by specifying the Y coordinate of the lower and upper edges of 
the spectra positions. The slitlet edges were automatically de- 
fined as the points where the intensity had dropped to 95% of 
the maximum level within the given slitlet screen flat spectrum. 
Both science and arc frames were cut. From this point on, the 
data reduction was identical to the reduction of longslit data. 

For each slitlet a 2D wavelength calibration was established 
and applied using the corresponding cut arc spectrum and stan- 
dard IRAF tasks. The rms of the fits was typically 0.05 A. 
Typically 30-40 arc lines well distributed over the wavelength 
range were used. 

Sky-subtraction was performed by first creating one- 
dimensional spatial profiles for each cut wavelength-calibrated 
science spectrum. Each profile was examined and suitable in- 
tervals along the spatial axis for which a fit to the sky level were 
decided. The sky was modelled as a constant at each pixel in 
the wavelength direction. 

One-dimensional spectra were extracted using the IRAF 
task APALL. The extraction aperture was chosen to be the 
FWHM of the spatial profile and within each aperture the flux 
was summed. 

The longslit spectra of the hot stars were reduced follow- 
ing a similar procedure to that of the MXU data. The issue of 
the two types of flat fields is iiTelevant here, since the extracted 
ID spectra were normalised to unity in the continuum. In the 
blue part of the spectrum, where these stars have intrinsic spec- 
tral features, the level was set to unity. The resulting telluric 
absorption correction spectrum was scaled and apphed to the 
galaxy spectra. 

As mentioned, slit spectra of the spectrophotometric stan- 
dard stars were taken at three positions: right, centre and left, 
covering three wavelength ranges: extreme blue, central and 
extreme red, respectively. The central spectrum has a substan- 
tial overlap with the two other spectra. However, when only the 



pixel-to-pixel flat field was applied, within the overlap regions 
the spectra had quite different shapes. We attributed this to the 
grism having a spectral response that depends on the position 
(angle) within the field of view. After the three spectra were di- 
vided by the screen flat field in which spectral shape was left in 
(cf. above), the spectra did match up and could be combined to 
form a single spectrum covering a wavelength range bracketing 
that for all galaxy spectra. From this spectrum the sensitivity 
function was derived and applied to the galaxy spectra that had 
been reduced using the same type of flat field. We thus use the 
fact that the position dependent spectral response of the grism 
is recorded both in the individual galaxy and star spectra and in 
the coiTesponding screen flat spectra, and when the former is 
divided by the latter, the spectral response of the grism cancels 
out, as long as the SED of the lamp stays constant. 

4. Cluster and galaxy redshifts 

4.1. Galaxy redshifts 

We present representative spectra for cluster members of the 
high-redshift cluster CI 1216 (figure and the intermediate- 
redshift cluster CI 1232 (figure^. We show the spectrum for 
the brightest cluster galaxy and spectra for an additional 5 pro- 
gressively fainter galaxies. 

Spectroscopic galaxy redshifts were measured using emis- 
sion lines where possible, in particular the [Oii]/l3727 line, or 
the most prominent absorption lines, e.g. Calcium K and H 
lines at 3934 A and 3968 A. Errors for galaxy redshifts were 
calculated using repeat galaxy observations acquired during 
spectroscopic runs 1 and 2. Our typical galaxy redshift er- 
ror was approximately 0.0003 (90 kms '), corresponding to 
50 and 60 kms"' rest frame at z = 0.8 and z - 0.5 re- 
spectively. We compared our redshift measurements with val- 
ue s determined u sing the Fourier Correlation Quotient method 
of lBendej(ll990l) and found a typical intrinsic uncertainty of 
0.00011 (18 kms"' atz = 0.8 rest frame; 22 kms ' atz = 0.5 
rest frame). 

Throughout this paper, a galaxy is considered to be a mem- 
ber if its redshift is within +(3xcr(,iuster) of the cluster redshift 
(cf. expression where crdustei is the cluster velocity disper- 
sion presented in Section|5] 

4.2. The spectroscopic catalogue 

We present the catalogue of spectroscopic observations of 5 
EDisCS clusters. The tables are available in electronic form at 
the CDS. The format of each table is illustrated in Table0] The 
tables provide the galaxy ID in column 1 ; object RA and Dec 
(J2000) in columns 2 and 3; I-band magnitude for an aperture 
of radius 1" in column 4; a spectroscopic redshift in column 
5; a membership flag in column 6; and a "targeting" flag in 
column 7. 

In column 1, a colon (:) indicates that multiple spectra were 
extracted for a single photometric object. Two redshifts have an 
associated ID with a colon: these two cases are shown in the ex- 
ample table (Table 1^. Our HST ACS image indicates that we 
are seeing two galaxies in projection. 3 galaxies were not in- 



o ^ PI 

o 
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EDCSNJiai6443-ia01429 
I 21.05 mag 




6000 6500 7000 

x(S) 



7500 




- EDCSNJ1216451-1 158493 
; I 21.16 mag 









- EDCSNJ121645S-1203134 
I 21.60 mag 



Hi; 



EDCSNJ1218434-120212E 
I 22.05 mag 



6000 



6500 



7000 

X(A) 



7500 



Figure 1. Example spectra (2 pixel binning) of galaxies in the cluster CI 1216 at z = 0.79. Galaxy IDs and 1" aperture I-band 
magnitudes are given for each spectrum. The bright emission line [OII]/l3727 (clearly visible in the right-hand panel figures) 
was used to determine the galaxy redshift where possible; otherwise the redshift was measured using the Calcium H and K lines 
(clearly visible in the interval 7000-7200A of the left-hand panel figures). All spectra are presented after correction for telluric 
absorption and before flux caUbration. 



Table 4. Example spectroscopic catalogue 



EDisCS object 


RA 


Dec 


I-band 


Spectroscopic 


Membership 


Targeting 








magnitude 


redshift 


flag 


flag 


EDCSNJ1054316-1147400 


10:54:31.62 


-11:47:40.0 


21.249 


0.6908 


1 


1 


EDCSNJ 1054254- 1147523: 


10:54:25.38 


-11:47:52.3 


21.380 


0.6977 


1 


1 


EDCSNJ1054254-1 147523: 


10:54:25.38 


-11:47:52.3 


21.380 


0.8385 





3 


EDCSNJ 1 054260- 1 148009 


10:54:26.02 


-11:48:00.9 


21.815 


0.6802 


3 


1 


EDCSXJ 1 054338- 1 145437 


10:54:33.81 


-11:45:43.7 


22.90 


0.6649 





1 


EDCSNJ 1054494- 1 24332 1 


10:54:49.36 


-12:43:32.1 


22.947 


0.6864:: 





1 


EDCSNJ1054504-1243410 


10:54:50.42 


-12:43:41.0 


21.471 


0.7309 


IB 


1 


EDCSNJ1054504-1243259 


10:54:50.44 


-12:43:25.9 


23.875 


0.225: 





3 


EDCSNJ1216431-1158113 


12:16:43.11 


-11:58:11.3 


22.503 


1.0579: 





1 


EDCSNJ1232389-1252170 


12:32:38.95 


-12:52:17.0 


21.523 


0.6809 





2 


EDCSNJ 1 232285- 1 254057 


12:32:28.54 


-12:54:05.7 


17.873 


0.0000 





4 


EDCSNJ1232322-1253003 


12:32:32.17 


-12:53:00.3 


20.824 


9.9999 





1 



eluded in the final versions of our photometric catalogues due 
to changes in object identification criteria. For these galaxies 
an empirical ID has been constructed using the galaxy RA and 
Dec coordinates, and a prefix of X instead of N; I-band magni- 
tudes are taken from the version of our photometric catalogues 
used to prepare the spectroscopic observations. 

In colunon 5, a colon indicates a redshift with a higher than 
average uncertainty, while a double colon (::) indicates a doubt- 
ful redshift. z = 0.0000 corresponds to stars, and z = 9.9999 is 
used where a redshift could not be determined. The uncertainty 
in the redshift estimate of the spectra listed with a colon was 



mostly due to a lower than usual S/N or sky features affecting 
the lines most suitable for the redshift measurements. This re- 
sulted in a less precise identification of the line centres. While 
the approximate redshift is confirmed by several lines in these 
cases, the exact redshift value cannot be established with the 
same accuracy as for the secure redshifts. Those cases with a 
double colon (6 spectra in total) are all low S/N spectra, where 
the redshift estimate had in general to be based on one single 
line in emission, whose redshift could not be confirmed by any 
other emission or absorption line in the spectrum. In all of these 
cases the reality of the observed line in emission was verified 
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Figure 2. Example spectra (2 pixel binning) of galaxies in the cluster CI 1232 at z = 0.54. Galaxy IDs and 1" aperture 1-band 
magnitudes are given for each spectrum. The bright emission line [011]/13727 (clearly visible in the two top right-hand panel 
figures) was used to determine the galaxy redshift where possible; otherwise the redshift was measured using the Ca H and K 
lines (clearly visible in the interval 6000-6200A of the left-hand panel figures). All spectra are presented after coiTection for 
telluric absorption and before flux calibration. 



by inspecting the 2D spectra. Only for one spectrum marked by 
a double colon was the redshift estimate based on the presence 
of a clear break in the continuum intensity, but no specific lines 
were identified without doubt due to the low S/N. 

The numeric membership flag in column 6 is "1" for cluster 
members (within Scrdustei of the cluster redshift); "2" and "3" 
for galaxies at 3crciu.stei to 4crciu.stei and AcTciuster to lOcr from the 
cluster redshift, respectively, that can be considered members 
of a supercluster region; and "0" in all other cases of fore- and 
back-ground non-members. The cluster CI 1054-12 has a sec- 
ondary structure lo-duster away in redshift from the main struc- 
ture (see Section 5 and Table 5). In this case a membership flag 
is Usted both for the main cluster CI 1054- 12A ("lA", "2 A", 
"3A", with meaning as above) and for the secondary structure 
CI 1054-12B ("IB", "2B", "3B"). Where a galaxy can have, 
for example, both a "3A" and "IB" flag, a "IB" flag has been 
assigned. 

Column 7 indicates whether the galaxy was a selected 
spectroscopic target, hence was targeted for spectroscopy as 
a possible member of the cluster (flag=l); if it was observed 
but believed to be a "field" (non-cluster) galaxy based on the 
photometric-redshift criteria (flag=2); if it is a galaxy that 
happened to fall into the slit together with a targeted galaxy 
(flag=3); and if it was targeted as a star to aid acquisition 
(flag=4). As described below, this information can be used to 
assess the success of our observing strategy and the effects of 
the selection criteria adopted. 



4.3. Success rate, completeness and potential 
selection biases 

It is important to establish a posteriori if the criteria used to 
select our spectroscopic targets have biased the final spectro- 
scopic cluster member dataset. As described in Section ITTl to 
maximise the number of targeted cluster members and avoid 
rejection of possible cluster galaxies, we adopted conservative 
criteria based on our photometric redshifts. To recap, our selec- 
tion criteria have rejected between 30% and 50% of the objects 
in the full photometric catalogues of our high redshift clusters, 
and approximately 35% of objects in the catalogue of the ~ 0.5 
cluster (CI 1232). 

The solid and dashed histograms in Figure |3 present the 
magnitude distribution of all objects detected by SExtractor in 
the EDisCS FORS2 field (i.e. our full photometric catalogues) 
and of potential spectroscopic targets decided after applying 
the selection criteria described in Section 12.11 The same plot 
illustrates the magnitude distribution of selected spectroscopic 
targets for which a spectrum was acquired (dotted histogram), 
and of selected spectroscopic targets which are cluster mem- 
bers according to their spectroscopic redshift (filled histogram). 

Figure |3] indicates that our membership success rate has 
been high. The fraction of objects targeted as cluster members 
and then confirmed to be members is approximately 75% in 
CI 1232, 60% in CI 1216, 50% in CI 1054-12 and CI 1054- 
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11, and 30% in CI 1040.^ The spectroscopic cluster catalogues 
typically include spectra for approximately 20% of all potential 
spectroscopic targets (~ 25% in CI 1216 and 35% in CI 1232). 
Redshift measurements were obtained for 91-98% of objects 
targeted as cluster members for a given set of cluster obser- 
vations. In cases where the spectrum did not yield a redshift 
the object was either of very low S/N spectra and among the 
faintest of the sample (generally fainter than / = 22.5), or rel- 
atively bright objects (/ ~ 20 or brighter) with a featureless 
spectrum that, judging from the shape of the red continuum, 
most probably are at a redshift much lower than our clusters. 

We recall that due to geometrical constraints in our mask 
design, 65 slits were assigned to objects rejected by our cluster 
member selection criteria. Importantly none of these objects 
was found to be a cluster member galaxy. Thus our adopted 
selection criteria have neither rejected real cluster members nor 
biased our spectroscopic catalogue of cluster members relative 
to a purely I-band selected sample (Section lzTl . 

5. Cluster redshifts and velocity dispersions 

In this section we present galaxy redshift histograms (Section 
15. U and cluster velocity dispersion measurements (Section]^) 
for our 5 EDisCS clusters. 

5.1. Redshift histograms 

In Figure0]we show redshift histograms obtained from the long 
spectroscopic observations of each cluster. The redshift distri- 
bution is clearly dominated by one main peak in redshift space 
at the cluster redshift. For CI 1054- 12, a secondary structure is 
evident very close to the main cluster peak. For CI 1232 there 
appears to be substructure in its velocity histogram. It is not im- 
mediately obvious that these structures are spatially offset from 
each other on the sky (see Section|6j. 

5.2. Cluster velocity dispersions 

We present the measurement of cluster velocity dispersion 
(o"ciiistei) for each cluster. 

A widely adopted estimator for c luster veloci t y disp ersions 
is the biweight statistic discussed bv Beer s et alJ (1 1990*) (here- 
inafter BFG90). For a typical dataset of 20 to 50 galaxies 
BFG90 show that the biweight scale estimator is a robust es- 
timator of distribution spread and demonstrate its insensitivity 
to outlying points. 

For each cluster, we estimate (Truster using biweight statis- 
tics: (Truster IS taken to be the biweight estimator of scale 
for an input redshift list and its asymmetric errorbars are es- 
timated from Monte Carlo bootstrap simulations (BFG90). 
Spectroscopic redshift measurements of uncertain quality (with 
a colon or double colon in Table 4) were rejected when estimat- 
ing these quantities. 

When computing the success rate, galaxies in CI 1054-12A and 
CI 1054-12B have been considered together, since the separation in 
redshift (0.0195) is far too small to be separated by any photometric- 
redshift technique. The lower success rate in CI 1040 is due to a higher 
contamination from back- and fore-ground galaxies. 



Cluster name 


^cluster 


'-^cluster — ^'-^clustei 


No. of 
members 1 


No. of 
members 11 


(1) 


(2) 


(3) 


(4) 


(5) 


CI 1040 


0.7043 


418 


30 


30 


CI 1054-11 


0.6972 




48 


49 


CI 1054-12 


0.7498 


504 :^^3 


35 


36 


CI 1216 


0.7943 


1018 


66 


67 


CI 1232 


0.5414 


1080 


52 


54 



Table 5. Cluster member numbers, redshifts and velocity dis- 
persions for our 5 EDisCS clusters. The cluster name is indi- 
cated in column 1. In column 2 we provide the cluster redshift. 
In column 3, cluster velocity dispersions (crdustei) are shown 
together with their errors in kms In column 4 we indicate 
the number of cluster members considered in the measurement 
of cluster velocity dispersion (see Section 15. 2> . In column 5 
the number of cluster members is revised to include members 
with redshifts of higher than average uncertainty and doubtful 
redshifts (see Sections f4.2l and l?!2t . 



An initial estimate of mean cluster redshift (zi„it) was made 
using redshifts judged to belong to the main peak of each clus- 
ter redshift histogram. 

Redshifts z,- within an interval of: 

(Zinn - 0.015) < Zi < (Zin„ + 0.015) (1) 

were used to obtain a first estimate of crdustei-- To transform to 
the cluster rest-frame, redshifts were multiplied by the speed of 
light, c, to create recessional velocities Vinit, ; the median of all 
these recessional velocities (Vmedian) was then found, and veloc- 
ities in the cluster reference frame were taken to be: 

/ ^median X 

Vresli - (Vini, - Vmedian)/(1 + )■ 

C 

An improved estimate of (Truster was taken to be the biweight 
estimator of scale for all these Vrest, - Its uncertainty was es- 
timated as the 1-cr bootstrap confidence interval. Next a 3-cr 
clipping was performed i.e. all input Vrest, outside the interval: 

— (3 X Cr cluster) ^ >^Ve.!r, < (3 X CT duster) (3) 

were rejected and the measurement of CTduster was repeated on 
the remaining velocity set. This step was iterated until expres- 
sion held true for all input Vrest, - The uncertainty of the final 
o"ciuster measurement was again estimated as the 1-cr bootstrap 
confidence interval. We define cluster members to be those 
galaxies for which expression Q holds after the final iteration. 

In Table |5] we give the results of this exercise. The cluster 
redshift given here, Zciuster, is the median redshift of all cluster 
members. 

6. Cluster substructure 

Figure |5] presents zoomed-in velocity histograms for our five 
clusters using 500 kms"' bins. In each panel, the red arrow 
marks the median recessional velocity of the cluster members 
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Figures. Magnitude distributions for our 5 clusters. The solid, largest histogram represents the magnitude distribution of all 
objects detected by SExtractor in the EDisCS fields (the full photometric catalogue in the area of sky sampled by the 2002 
spectroscopy). The dashed histogram is the distribution of objects that were retained as potential spectroscopic targets within 
the area of sky sampled by the 2002 spectroscopy after applying the selection criteria described in Section ITTI The magnitude 
distribution of selected spectroscopic targets for which a spectrum was obtained is shown as a dotted histogram, while that of all 
selected spectroscopic targets confirmed to be cluster members is shown as a filled histogram. 
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Figure 4. Redshift histograms based on the long spectroscopic exposures of our 5 EDisCS clusters. An inset panel in each plot 
shows a zoomed histogram for the redshift interval Zciuster - 0.04 < z < Zciuster + 0.04, where the cluster redshift Zciuster is indicated 
by a dotted line. Dashed lines indicate the interval Zciuster - 0.015 < z < Zciuster + 0.015. Redshift bins are 0.00625 (1874 kms"' 
observed frame) for the main plot and 0.000625 for the inset plot. In all cases the redshift distribution is dominated by one peak 
at the cluster redshift. 
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(i.e. the cluster redshift defined above), and the blue arrow 
marks the recessional velocity of the BCG. A Gaussian centred 
at the median recession velocity and of width corresponding 
to the cluster velocity dispersion is plotted as reference. The 
velocity histograms show different levels of deviation from a 
Gaussian, providing evidence of complex dynamical structure 
and suggesting that some clusters are far from equilibrium. In 
this Section, we assess substructure within each cluster first 
based on XY position maps of cluster members for different 
redshift intervals (Section 16. 1> and then using the Dressler- 
Shectman test (Section l6^ . 

6.1. XY position diagrams 

In Figure |6l we present XY position maps for galaxies within 
the I-band magnitude range for which our spectroscopy was 
completed and we consider as possible cluster members on 
the basis of their photo-z's. Objects identified as non-members 
from our spectroscopy are excluded. For each cluster the posi- 
tion of the BCG is indicated by both a black cross and a colour 
symbol. We study cluster structure in three dimensions by sub- 
dividing cluster members into three groups corresponding to 
different redshift intervals. Members indicated by green circles 
are all galaxies with redshift z, in the interval 

(2clu.ster ~ Zo-) ^ Zi < (Zcluster + Zcr) 

where Za- = S^ish. x (1 -i- z) is the redshift interval correspond- 
ing to the cluster velocity dispersion (Truster presented for each 
cluster in Section ls!2l 

Similarly blue triangles indicate cluster members with red- 
shift Zi such that 

Zi < (Zcluster ~ Zcr) 

and red squares, cluster members with redshift z, such that 

Zi > (^cluster + Zcr)- 

For CI 1054-12, some clustering of galaxies indicated by 
red and green symbols may be present close to the cluster cen- 
tre. For CI 1216, there are two clear substructures, one given by 
galaxies with red and green symbols at cluster centre and an- 
other less concentrated structure offset towards low X and high 
Y (northeast of the BCG), shown mainly by blue symbols. For 
CI 1232 an extended substructure indicated mainly by red sym- 
bols is clear towards the cluster centre. No clear spatial segre- 
gation is apparent for any redshift intervals in CI 1054-1 1 and 
CI 1040. 

6.2. The Dressier-Shectman test 

More robust evidence for substru cture requires its significance 
and its properties to be quantified. lDressler & Shectmanlj 19881) 
devised a statistical test that uses galaxy velocities and posi- 
tions to constrain cluster substructure in three dimensions. The 
principle of the method is the following. Starting from a list 
of cluster members with measured positions and velocities, for 
each galaxy one finds the ten listed nearest neighbours on the 
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Cluster name 




A 


P 


CI 1040 


30 


55.11 


0.08 


C11054-11 


48 


71.77 


0.64 


CI 1054-12 


37 


108.87 


0.14 


C11216 


67 


135.58 


0.05 


CI 1232 


54 


96.57 


0.01 



Table 6. Number of spectroscopic members (A^;,), value of the 
A statistic, and the measure of substructuring significance P, 
for each of the clusters used in the present analysis (see text for 
details). 



sky. The local mean velocity and velocity dispersion are com- 
puted from this sample of 1 1 galaxies. These quantities, defined 
for each galaxy in the list, are then compared to the global clus- 
ter mean and velocity dispersion using the parameter 6 defined 
to be: 

6^ = [(Viocal - vf + (o-iocal - 0")^] (4) 

where cr and v are the cluster velocity dispersion and the clus- 
ter redshift respectively, as given for each cluster in Section|5l 
while cTjocai and viocai are the quantities defined locally for each 
gal axy. 

IPressler & S hectmanI (Il988l) define the cumulative devia- 
tion A as the sum of the 6 for all the cluster members A^g. Note 
that the A statistic is similar to a;^'^: if the cluster velocity dis- 
tribution is close to Gaussian and local variations are just sam- 
pling fluctuations, then A will be of the order of A^g. 

This technique does not allow a direct identification of 
galaxies belonging to a detected substructure. The position of 
substructures can however be identified by plotting the distri- 
bution of the galaxies on the sky using symbols whose size is 
proportional to the parameter 6, thus quantifying the local de- 
viation from the global kinematics of the cluster. Such an anal- 
ysis is shown for our 5 clusters in Figure In each panel, the 
size of the symbols is proportional to e^^^ and a cross symbol 
marks the position of the BCG. The circle symbols are coloured 
according to the criteria described in Section lOl To first or- 
der the presence of many large circles in a given area should 
correspond to substructure, i.e. a correlation both spatially and 
kinematically of separate structures. Our plots in Figure in- 
dicate that complex dynamical structure may be present for all 
clusters although significant substructuring is detected for only 
CI 1216 and CI 1232. 

The A statistic can be used to give a quantitative estimate of 
the significance of the putative substructure. We have computed 
a set of 1000 Monte Carlo simulations for each cluster by ran- 
domly re-shuffling the velocities of the cluster members. This 
removes any significant 3D substructure. The statistical signif- 
icance of the o bserved substructure can th en be quantified by 
the fraction P dPress ler & Shectman' '1988^ of the simulations 
that yield A values larger than the observed one where a small 
value of P corresponds to a high significance. In Table|6lwe list 
the number of spectroscopic members, the measured A statistic, 
and the measure of substructuring significance P. 
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Figure 5. Velocity histograms for the five clusters discussed in this paper The red open symbol arrows mark the median reces- 
sional velocity of all the spectroscopically-confirmed cluster members, while the blue filled arrows mark the recessional velocity 
of the BCGs. The red Gaussians correspond to the measured velocity dispersion. 



For two of our clusters (CI 1232 and CI 1216), substruc- 
ture is detected by the DS test with more than 95% confidence. 
This agrees with the spatial segregation of galaxies in differ- 
ent velocity intervals seen in Figure|6land with the existence of 
peaks redward (CI 1232) or blueward (CI 1216) of the cluster 
redshift in the velocity histograms (Figure|5}- In CI 1232 the 
red-shifted galaxies are mainly projected near the cluster cen- 
tre, while in CI 1216 the blue-shifted galaxies are offset North- 
East of the BCG. In two clusters (CI 1040 and CI 1054-12), 
substructure appears to be present, based on the DS probabil- 
ity, but the available data do not yet provide firm statistical 
evidence. No evidence for substructure is found in CI 1054- 



1 1 . Interestingly, the two clusters with the clearest substructure 
are the ones with the largest velocity dispersions, both exceed- 
ing 1000 kms At these redshifts, equilibrium clusters with 
such high velocity dispersions are rare according to the stan- 
dard hierarchical ACDM scenario (e.g. Eke et al1 ll996l) . The 
substructure analysis suggests that both clusters have consider- 
able depth in redshift, indicating that they may still be collaps- 
ing along or close to the line-of-sight. One should therefore be 
extremely cautious when using their measured velocity disper- 
sions to estimate their masses. Such non-virialised structures 
are not unexpected in ACDM models cEvrard et al.-2002) . 
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Figure 6. Plots of galaxy positions. Axis scales are in units of CCD pixels (0.2"/pixel). The cluster name is indicated at the top 
of each plot. Large open symbols show the positions of cluster members: green circles indicate galaxies within + one Zcr of the 
cluster redshift, red squares indicate galaxies whose redshift is more than one Za- higher than Zduster, while blue triangles indicate 
galaxies whose redshift is lower than Zdustei by more than one Zcr where Za- is the redshift interval corresponding to the cluster 
velocity dispersion in observed frame. Small circles are galaxies without spectroscopy selected from our 1-band photometric 
catalogues to be possible cluster members. The position of the brightest cluster galaxy (BCG) is indicated by a black cross 
symbol. A physical distance of 1 Mpc (assuming Q = 0.3, A = 0.7 and h - 0.7) is shown in each plot. 
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Figure 7. Dressler-Shectman (DS) analysis. The radii of the plotted circles is equal to where S is the DS measurement of local 
deviation from the global velocity dispersion and mean recession velocity. Green circles marked by a solid line indicate galaxies 
with redshift Zi such that {zduster - Za-) < Zi < (zduster + Zo-), red circles indicated by hashed areas, galaxies with Zi > (zdusier + Za-), 
and blue circles indicated by a dotted line, z, < {zduster - o'duster)- The measure of substructuring significance P is provided at the 
top right hand corner of plots for each cluster. 
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7. Summary 

We present spectroscopic observations of 5 clusters, 4 at red- 
shift of 0.7-0.8, 1 at redshift ~ 0.5, completed as part of the 
ESO Distant Cluster Survey (EDisCS). We demonstrate how 
our target selection strategy has maximised our cluster member 
dataset. We present spectroscopic catalogues of positions, red- 
shifts and 1-band magnitudes for galaxies observed in the fields 
of our clusters. We have acquired data for 236 cluster members 
and the number of cluster members per cluster ranges between 
30 and 67. Our spectroscopic success rate at targeting mem- 
bers was 30-75%, averaging 50%. Our cluster velocity disper- 
sions measured using a robust biweight estimator range from 
400 to 1 100 km s ' . We find significant evidence for substruc- 
ture in two clusters, one at z ~ 0.79 and the other at z ~ 0.54. 
Both clusters have velocity dispersions exceeding 1000 kms"' 
and are unlikely to be fully virialised: their velocity dispersions 
may therefore not provide a robust estimate of their mass. The 
properties of our 5 EDisCS clusters cover a range of veloc- 
ity dispersions, richness and substructuring typical of our sur- 
vey. Our full spectroscopic dataset will allow a comprehensive 
study of galaxy evolution as a function of redshift and cluster 
environment. 
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Appendix A: Cosmic ray removal algorithm 

For the 'long masks' typically 4 or 8 exposures of 
30 minutes were availa ble and the following method was used 
(lMilvang-J ensen'2003': Milvang-Jensen & Aragon-Salamanca, 
in prep.). Consider the example of 8 exposures. To be able to 
flag deviating values among the 8 values available at a given 
pixel position we must calculate the expected standard devia- 
tion (hereinafter 'sigma') in the absence of cosmics. One es- 
timate of sigma can be computed from the median using the 
CCD noise model (photon noise and read-out noise). Since 
the numerous sky lines have large intrinsic intensity variations 
this sigma can however significantly underestimate the actual 
frame-to-frame variation. The sigma must therefore be mea- 
sured from the data. The number of pixels alfected by cosmics 
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in a given 30 minute exposure is approximately 0.2%. It fol- 
lows that for a given pixel position at most 2 values are ex- 
pected to be affected by cosmics, and thus the lowest 6 values 
can be used to compute a sigma. As a statistic we used the RMS 
(run 2) or the MAD (mean absolute deviation from the sample 
mean) (run 3); the results were found to be similar. This sigma 
estimate is biased low with respect to the intrinsic value (still 
in the absence of cosmics). The exact factor was found from 
simulations using a normal parent distribution. The estimates 
were then multiplied by this factor to give a 'measured' sigma 
estimate, e.g. cTmeas. = 1.79-MAD(lowest6of8). Thisestimate 
incorporates the large variations due to the varying strength of 
the sky lines, while giving a value close to that predicted by 
the CCD noise model in regions of the spectra not affected by 
sky lines. However, due to small number statistics (~ V6 un- 
certainty) the measured sigma can occasionally be too small. 
Therefore, the finiil sigma was taken to be the maximum of 
the measured sigma and the CCD noise model sigma. Cosmics 
in the individual frames could now be identified as values ex- 
ceeding the median image by more than A'^ sigma. The factor 
was chosen manually on a frame by frame basis in such a 
way that very few good (i.e. non-cosmic) pixels were flagged. 
Usually A'^ = 4-5 worked. If however a series of 8 (say) mask 
exposures was started at the beginning of the night, sky lines 
in the first exposures could be much brighter than in the re- 
mainder in a way not predictable using the measured sigma: 
in this case a larger threshold was required. For run 2 the 
cosmic maps were used to clean the individual exposures us- 
ing Unear interpolation within each frame. Following this the 
exposures were combined (i.e. averaged). For run 3 the cos- 
mic masks were used to exclude cosmics when combining the 
frames. It should be noted that the stability of F0RS2 is ex- 
cellent and the shift in wavelength and spatial direction during 
4 hours (8 exposures) of uninterrupted observations is small 
(< 0.5 pixel), allowing frames to be combined without shifts. 



